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Abstract. A novel quantum dynamical model based on the dissipative quantum dynamics of open 
quantum systems is presented. It allows the treatment of both deep-inelastic processes and quantum 
tunneling (fusion) within a fully quantum mechanical coupled-channels approach. Model calcu- 
lations show the transition from pure state (coherent) to mixed state (decoherent and dissipative) 
dynamics during a near-barrier nuclear collision. Energy dissipation, due to irreversible decay of 
giant-dipole excitations of the interacting nuclei, results in hindrance of quantum tunneling. 
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INTRODUCTION 

Stationary state coupled-channels approaches have been very successful fP] in explain- 
ing several collision observables. However, there are still unsolved problems. For in- 
stance, the inability to describe elastic scattering and fusion measurements simultane- 
ously 10, 111] and, related, the more recent failure to describe in a physically consistent 
way the below-barrier quantum tunneling and above-barrier fusion yields [4]. 

These problems may be caused by the neglect of important physical processes (e.g., 
deep-inelastic) which cannot be treated within (standard) coupled-channels models. 
Measurements have shown that deep-inelastic processes occur even at sub-barrier in- 
cident energies [jSt], in competition with the process of quantum tunneling, and thus 
fusion. Energy dissipation associated with the deep-inelastic mechanism could thus play 
a significant role in the inhibition of tunneling at deep sub-barrier energies. This can 
also change the yield of direct reaction processes, including elastic and quasi-elastic 
channels. 

The understanding of this complex interplay, at near- and below barrier energies, re- 
quires a dynamical model which can describe coupling assisted tunneling with dissi- 
pation. Neither existing models of fusion nor of deep-inelastic scattering can address 
both energy dissipation and quantum tunneling. Quantum mechanical coupled-channels 
models describe tunneling without energy dissipation whilst approaches to direct 

damped collisions treat the relative motion of the nuclei classically [TJ- 



We here report on a novel coupled-channels density matrix approach fs*] that over- 
comes these difficulties. We exploit the Lindblad axiomatic approach [9, 10] for open 
quantum systems, which is based on the concept of quantum semigroups and completely 
positive mappings. We refer to Ref. [8] for a discussion of the suitability of the Lind- 
blad theory for the treatment of low-energy collision dynamics. The coupled-channels 
description is formulated with Lindblad's equation for a reduced density matrix [9,. .1QJ . 
It describes the dynamical evolution of the reduced system (comprising the relative mo- 
tion of the nuclei plus selected, intrinsic collective excitations) that irreversibly interacts 
with two (model) "environments". Firstly, an environment inside the Coulomb barrier, 
which is related to the complexity of compound nucleus states. Secondly, one with a 
long range, associated with decay out of short lived (compared to the reaction time) in- 
ternal vibrational states, e.g. the giant dipole resonance (GDR) of the colliding nuclei. 
Model calculations show that damping of the GDR results in quantum decoherence and 
energy loss as the nuclei overlap, inhibiting tunneling, and thus fusion. 

COUPLED CHANNELS DENSITY MATRIX APPROACH 

The model calculations are carried out using the scenario presented in Figure [IJ The 
basis will comprise two asymptotic states (coupled-channels) 1 1) and |2). Channel 1 1) is 
the (ground states) entrance channel and is coupled to an inelastic state |2) by a coupling 
interaction V12. Two distinct sources of irreversibility are also considered, modelled by 
two auxiliary (environment) states \X) and \Y). The first environmental coupling de- 
scribes capture by the potential pocket inside the fusion barrier. This simulates the irre- 
versible and dissipative excitations associated with the evolution from the two separate 
nuclei to a compound nuclear system. In a stationary states approach this loss of flux 
is approximated by imposing an imaginary potential —iW{r),W{r) > 0, or an ingoing- 
wave boundary condition at distances well inside the barrier. Here, these transitions are 
described by an auxiliary state \X), to which all other states \j) couple, modelled llllll 
by a Lindblad operator ^xj = VY^\^){j\ - The absorption rate to state \X) is given by 
Y'^ = W{r)/h where W{r) is taken as a Fermi function with depth 10 MeV and dif- 
fuseness 0.1 fm, located at the pocket radius of the nucleus-nucleus potential, ^ 7 fm. 
This choice guarantees complete absorption inside the pocket. The fusion probability is 
defined as the probability accumulating in this state \X). 

The second environment, whose explicit treatment will be seen to be the most signif- 
icant at lower energies, is associated with the irreversible decay out of intrinsic excita- 
tions of the colliding nuclei. Such decays are independent of the dynamical couplings. 
Specifically, we will associate the only excited coupled channel state |2) with the GDR 
excitation. We then introduce a second auxiliary state \ Y), representing the bath of states 
in which the GDR is embedded, and to which only the GDR excitation |2) is coupled. 

Thus, \Y) and/or \X) supplement the two intrinsic states |1) and |2) that comprise the 
two coupled channels. Both of the auxiliary states refer to complex excitation modes of 
the nuclei, associated with nucleonic degrees of freedom and compound nucleus states, 
respectively. They provide intuitive and formal channels [11] for describing irreversible 
coupling and loss of probability from the system to these environments, couplings that 
enter only through Lindblad's dissipative Liouvillian [8] . \Y) is also assumed to couple to 
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|X) - Environment associated with compound nucleus states 
Cxj = VT^I^)(j| where 7'^'^ = W{r)/h 
\Y) - Environment associated with nucleonic degrees of freedom 
Cyj = \/TYj\Y) (j\ where Tyj is a constant 

FIGURE 1. (Color online) Schematic picture of the model scenario, and '^yj are physically 
motivated Lindblad's operators for the two dissipative channels \X) and \Y), respectively. See text for 
further details. 



\X) at the appropriate range of separations. Probability accumulating in state \Y) outside 
of this \X) pocket may be identified with deep-inelastic processes. 

Dynamical calculations proceed as follows. The reduced density matrix results from 
the representation of the reduced density operator p (?) in an asymptotic (product) ba- 
sis of coordinate states and intrinsic energy states of the individual nuclei. Its dynamical 
evolution is dictated by a set of Lindblad's coupled equations [8] for the matrix elements, 
which defines a time-dependent coupled-channels problem with an initial value. The 
initial pure-state density matrix, describing the well- separated nuclei (in their ground 
states) with a wave-packet in their relative motion, has Tr[p^] = 1. The purity of this 
state, conserved under unitary time evolution, is destroyed (Tr[p^] < 1) if the environ- 
ment causes a loss of quantum coherence. The degree of decoherence is quantified by 
a loss of purity of the density matrix. Having solved the time-evolution of the density- 
matrix numerically [8], the expectation value of reaction observables ^ is given by eval- 
uating the appropriate trace {S'{t)) = Tr[^p(f)]. In this way, probabilities of reaction 
channels can be determined. 

We have performed test calculations for a head-on collision of ^^O -1- ^^^Sm at energies 



below the nominal Coulomb barrier (Vb = 61.1 MeV). The Hamiltonian of the reduced 
system is the same as that in the fusion model implemented in CCFULL [6]. The form 
of the bare nuclear potential between the two nuclei, consistent with the stated Vb, is a 
Woods-Saxon potential with (Vq, ro, ao) = (— 105. 1 MeV, 1 . 1 fm, 0.75 fm). The Coulomb 
potential was that for two point charges. The ^^O projectile was taken to be inert and 
the ^^^Sm target was allowed to be excited to a GDR vibrational state. The dynamical 
nuclear coupling of the ground state 1 1) to the vibrational state |2), with excitation energy 
El- = 15 MeV, has a macroscopic deformed Woods-Saxon form with a deformation 
parameter of /3i = 0.2. 

The time step for the density-matrix propagation was At = 10^^^ s, and the radial 
grid (r = — 250 fm) was evenly spaced with M = 512 points. The relative motion 
of the two nuclei in the entrance channel 1 1 ) was described by a minimal-uncertainty 
Gaussian wave-packet, with width Oq = 20 fm, initially centered at r = 150 fm, and was 
boosted towards the target with the appropriate average kinetic energy for the entrance 
channel energy Eq required. The FWHM energy spread of the wave-packet is ~ 3%. 
The numerical accuracy of the time evolution was checked using a fully coherent, time- 
dependent calculation, excluding coupling to states \X) and \ Y).lt was confirmed that the 
normalisation and purity of the density-matrix, Tr[p] = Tr[p^] = 1, and the expectation 
value of the system energy Tr[^p] were maintained with high accuracy over the required 
number of time steps, typically 700 for the full duration of the collision. 

The importance of the two, spatially distinct, sources of environment couplings were 
studied. Calculations were first performed only including the effect of coupling of the 
intrinsic coupled channels 1 1) and |2) to the capture state \X). Calculations were carried 
out for Eq = 45, 50, 55 and 60 MeV incident energy. The calculated state purity Tr[p^] 
and the energy dissipation Tr[^(po — p)] post the collision (after 700 time steps) are 
shown in the left panels in Table [TJ For sufficiently sub-barrier energies, £0^55 MeV, it 
is evident that time-evolution in the presence of state \X) essentially maintains coherence 
and is non-dissipative. There is however loss of purity and dissipation at the highest 
energy. It is interesting therefore to compare the density-matrix and the predictions of 
CCFULL (that uses stationary Schrodinger's dynamics with an ingoing wave boundary 
condition). This is done here only for calculations of the tunneling probability P{Eq) , in a 
relative s-wave, shown in Figure[2l These comparisons, of necessity, require convolution 
of the £ — partial wave penetrabilities To{E) from CCFULL with the energy distribution 
/{E.Eq) of the chosen initial wave packet. That is, P{Eq) = J dE f{E,Eo) Tq{E). The 
P{Eo), shown as a function of Eq/Vb in Figure l2l are in very good agreement showing 
the appropriateness of stationary state coupled-channels calculations for this observable 
within the dynamical scheme of states |2) and \X). It is our contention that the 
dissipation associated with state \X), while significant at 60 MeV, is strongly localised 
inside the barrier and thus does not impact upon the barrier penetrability. We will now 
show that the same is not true for the more spatially-extended dissipation due to the 
GDR decay environment \Y). 

The treatment of the irreversible GDR decay (with a spreading width of 6 MeV) to the 
bath of surrounding complex states (represented by \Y)) was included by switching on 
the coupling of the intrinsic inelastic state |2) to \Y). Unlike the coupling to \X), a major 
part of the inelastic excitation of the system gives access to \Y) before the wave packet 
encounters the fusion barrier. The onset of decoherence, the purity of the density matrix. 



TABLE 1. The calculated density matrix purity Tr[p^] and energy 
loss AE —Ti[H{j5o — p )] following time-evolution (for 700 time steps) 
when including only the state \X) (left entries) and both states \X) 
and |F) (right entries) environmental couplings. The GDR coupling 
strength used was j3i = 0.2. 



State \X) I States \X) and \Y) 



Eo (MeV) 


Tr[p2] 


AE (MeV) 


Tr[p2] 


AE (MeV) 


45 


1.0000 


0.0004 


0.9196 


1.8718 


50 


1.0000 


0.0004 


0.8977 


2.6744 


55 


0.9996 


0.0109 


0.8759 


3.6100 


60 


0.6067 


14.862 


0.5127 


18.908 
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FIGURE 2. (Color online) The energy dependence of the s-wave tunneling probability calculated with 
the density matrix (solid points) and the coupled-channels CCFULL methods (full line). See text for further 
details. 

and the associated energy dissipation are shown in the right hand entries in Table [IJ 

The probability trapped under the fusion barrier is associated with GDR collective vi- 
brational energy being irreversibly removed from the coherent dynamics into surround- 
ing bath states (heat). This is then no longer available for relative motion, or tunneling. 
Such energy loss can be correlated with deep inelastic processes, seen experimentally, 
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FIGURE 3. (Color online) Time-dependence of the probability trapped in \X) for £0 = 45 MeV. The 
full curve includes states |2) and \X). The dotted curve adds the irreversible decay of |2) to \Y). The 
calculations are for )3i = 0.2. 

that compete with fusion in reactions involving heavy nuclei fs*]. 

Figure [3] shows the time evolution of the probability trapped in the potential pocket, 
state \X), for Eq = 45 MeV. We comment that, when including the inelastic channel 
|2) but not \Y), the nucleus-nucleus potential renormalization leads to the expected 
enhanced penetrability from the inelastic channel coupling, compared to the purely 
elastic (|1) plus \X)) calculation. The decoherent dynamics due only to environment 
\X) gives the (full curve). By comparison, the calculation that also includes the GDR 
doorway-state decay to \Y) leads to a suppression (dotted curve and arrow) of the 
population of state \X). Additional irreversible processes other than excitation of the 
GDR are also likely to contribute to the deep-inelastic yield, such as complicated multi- 
nucleon transfers [12]. To simulate these very simply, the assumed state |1) to |2) 
coupling strength was increased. Figure |4] shows the dependence of the calculated 
tunneling suppression on the assumed j8i strength for Eo = 45, where we note that larger 
/3i result in both an increase in the strength and the range of the coupling formfactor to 
the inelastic state |2). 
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FIGURE 4. Calculated suppression of the probability trapped in \X) as a function of the assumed j3i 
value for Eq = 45 MeV. 



SUMMARY 

We have reported on a quantum dynamical approach based on time-propagation of a 
coupled-channels density matrix. Both deep-inelastic processes and quantum tunneling 
(fusion) can be treated within this fully quantal framework. It describes the transition 
from pure state (coherent) to mixed state (decoherent and dissipative) dynamics during 
a nuclear collision. The development provides a significant step towards an improved 
theoretical understanding of low-energy collision dynamics, as the calculations exhibit 
both quantum decoherence and energy dissipation. Effects of decoherence and dissi- 
pation on collision dynamics can be manifested at distances outside the fusion barrier 
radius, resulting in suppression of the quantum tunneling probability. This may have 
major implications for understanding current problems in near-barrier reaction dynam- 
ics [0, [3, 4], including the sub-barrier fusion hindrance phenomenon |13]. More com- 
plete calculations and detailed consideration of other processes, such as multi-nucleon 
or cluster transfer reactions, are required to confront measurements. 
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